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Abstract 

The Majorana Experiment will use arrays of enriched HPGe detectors to search for the neutrinoless double-beta decay of 76 Ge. 
Such a decay, if found, would show lepton-number violation and confirm the Majorana nature of the neutrino. Searches for such 
rare events are hindered by obscuring backgrounds which must be understood and mitigated as much as possible. A potentially 
important background contribution to this and other double-beta decay experiments could come from decays of alpha-emitting 
isotopes in the 232 Th and 238 U decay chains on or near the surfaces of the detectors. An alpha particle emitted external to an HPGe 
crystal can lose energy before entering the active region of the detector, either in some external-bulk material or within the dead 
region of the crystal. The measured energy of the event will only correspond to a partial amount of the total kinetic energy of the 
alpha and might obscure the signal from neutrinoless double-beta decay. A test stand was built and measurements were performed 
to quantitatively assess this background. We present results from these measurements and compare them to simulations using 
Geant4. These results are then used to measure the alpha backgrounds in an underground detector in situ. We also make estimates 
of surface contamination tolerances for double-beta decay experiments using solid-state detectors. 



1. Introduction 

Observation of neutrinoless double-beta decay (Ov/3/3) rep- 
resents the best chance for discovering the nature of the neu- 
trino (Majorana or Dirac) HI Q El B] ■ Experimental searches 
using high-purity germanium (HPGe) detectors — enriched to 
86% 76 Ge — have demonstrated the most stringent half-life lim- 
its for this decay (T°J 2 > 1.9 x 10 25 y [5]) to date. Future ex- 
periments with tonne-years of exposure should be able to reach 
, ^ half-life limits greater than 10 27 y, but only through concerted 
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efforts in understanding and reducing backgrounds. Current ex- 
periments using 76 Ge, such as Majorana||6[ [H and GERDA 
Phase IlHm, 

will attempt to reach background levels of ~ 
4 background events/tonne-year in the Ov/3/3 region-of-interest 
around 2039 keV to demonstrate the viability of future tonne- 
scale experiment^ Achieving these background goals will re- 
quire ultra-clean materials, sufficient shielding, sophisticated 
background-rejection techniques, and a deep-underground set- 
ting. 
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3 A background of 4 events/tonne-year in the Majorana and GERDA exper- 
iments will scale to ~ 1 event/tonne-year in a tonne-sized experiment. 



The signature of Ov/3/3 is a peak at the Q-value of the de- 
cay, corresponding to the sum of the kinetic energies of the two 
emitted electrons. Alpha particles, emitted from decays in the 
232 Th and 238 U natural decay chains with energies of 3.9-8.8 
MeV, can lose kinetic energy before entering the active region 
of an HPGe crystal. These degraded alphas can result in a con- 
tinuum of events, obscuring a possible signal from Ov/3/3. 

HPGe diode detectors have n + and p + ohmic contacts for 
bias and signal connection. The n + layer, created with diffused 
lithium ions, is typically 0.5 - 1 mm in depth. The p + layer 
is created by implanting boron ions, resulting in a dead layer 
on the order of ~ 0.5pm. These form regions within the de- 
tector that are insensitive to ionizing radiation, where energy 
lost by a particle will not be registered. Because alphas can- 
not traverse the thick n + layer, a detector's susceptibility to al- 
phas is dependent upon the amount of surface with a p + con- 
tact. HPGe detectors are characterized by their dopants, be- 
ing either n-type or p-type. The dopants then determine the 
type of ohmic contacts for the crystal; n-type detectors have n + 
central contacts while p-types have p + . The outer surface con- 
tacts are reversed for each detector type. The diagram in Fig. 
Q] shows a depiction of these detector types with exaggerated 
thicknesses. The thin outer contact means that n-type detec- 
tors have more alpha-susceptible surface area than p-types. A 
p-type with a point contact (p-pc, 111111 ) will be used in the Ma- 
jorana Demonstrator||7| and Phase II of GERDA 1 12 ] . The re- 
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gion between the n + and p + layers is typically passivated and is 
characterized by incomplete and hard to predict charge collec- 
tion (shown as dashed in the figure). The response of a detector 
to alphas incident on this surface is still an unknown. The p- 
pc detectors used in the Majorana Demonstrator will have a 
minimal amount of this type of surface. 



n-Type 



p-Type 



p-pc 



Figure 1: Cross-section diagrams of different HPGe detectors. 
The thin outline represents the p + (thin) dead layer and the thick 
outline represents the n + dead layer. The n- and p-type detectors 
are of the common semi-coax design. The p-type point contact 
detector has far less thin surface than either the p- or n-types. 
The area on the bottom of the crystals (dashed line) represents 
the passivated surface that insulates the n + and p + contacts. 

Broadly speaking, alpha backgrounds for QvBB experiments 
with HPGe detectors fall into two categories, depending on 
where the degradation in kinetic energy occurs: surface alphas 
and external-bulk alphas. 

Surface-type alpha events are characterized by decays at the 
surface of the crystal, with the alpha losing kinetic energy 
within the dead region of the crystal. Only alphas that enter 
the crystal at extremely shallow incidence angles (with respect 
to the normal of the crystal surface) will lose an appreciable 
amount of energy within this thin dead region; alphas entering 
the crystal at or near normal to the surface will lose a negligible 
amount of energy before entering the active region of the crys- 
tal, resulting in minimal energy loss and a peak structure when 
measured by the detector. The classic example of this back- 
ground arises from exposure to 222 Rn. The decay of 222 Rn (3.8 
days) and its subsequent daughters 2I8 Po, 214 Pb, 214 Bi, 214 Po 
eventually lead to 2I() Pb. Along the way, these daughters of 
222 Rn can implant onto surfaces llOl . 210 Pb decays to 21() Bi, and 
then 210 Po, which emits a 5.3 MeV alpha upon its decay. The 
relatively-long half-life of 210 Pb (22 years) means that even a 
brief exposure of a detector or its surroundings can lead to a 
steady supply of 5.3 MeV alpha backgrounds ! l4 1 1 311 . 

The other category of alpha backgrounds, external-bulk al- 
phas, originate in bulk materials external to the crystal, e.g. 
a contact pin or detector mount. Energy loss of the alpha oc- 
curs in this bulk material before it hits the detector surface and 
depends on the amount of external material the alpha travels 
through. The result is a broad continuum of events with no 
peak structure. 

This paper presents work to quantify the role of alpha back- 
grounds on Qv/3/3 experiments using HPGe detectors, in Sec. 
[2] a test stand that was designed and built to study alpha back- 
grounds is described. Also in this section a model of the re- 



sponse of an HPGe detector to alpha decays is discussed and 
applied to the test stand data. In Sec. [3]this model is applied to 
data taken with an in situ detector located underground. These 
results are discussed in the context of backgrounds to Ov/3/3 ex- 
periments using HPGe detectors in Sec. |4] 

2. Background Test Stand 

An n-type, semi-coax HPGe detector manufactured by 
ORTEC[15] was modified to facilitate measurement of an al- 
pha source on the surface of an HPGe crystal. The original de- 
tector was enclosed within an aluminum outer vacuum cryostat 
that sealed against the base of the detector canister. The base 
is mounted securely to the cold finger and contains all electri- 
cal feedthroughs. The inside of the canister contains the upper 
cold finger and cold plate, front-end electronics for the pream- 
plifier, and the crystal. The crystal mount consists of a smaller 
aluminum cup attached to the cold finger with the HPGe crys- 
tal affixed inside this mounting cup. A thin sheet of aluminized 
mylar that originally covered the crystal's front face was re- 
placed with a crystal cap fabricated from aluminum. The exter- 
nal vacuum enclosure — sealed against the detector base with an 
O-ring — was replaced with a larger, custom cryostat enclosure 
for ease of source transfer and for external-source manipulation. 
The test stand was operated in two separate modes to study both 
surface-type and bulk-type alpha backgrounds. 




Figure 2: Side-view of top portion of modified detector cryostat 
enclosure. A rotational feedthrough (a) is built into the the mod- 
ified cryostat enclosure (b). The feedthrough rotates a source 
(c) in a circle above the collimation plate (d). The collimation 
plate (Fig. [3]) rests on the crystal mounting cup above the HPGe 
crystal (e) and provides a collimated source of alphas from the 
source to illuminate the crystal face. 



2.1. Surface Background Model 

For studying surface-type backgrounds, the test stand (see 
schematic, Fig. |2]i employed a rotational feedthrough (a) 
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Table 1: The alpha decay of 241 Am has five prominent alphas 
and two prominent gamma lines from the decay to 237 Np (22.35 
keV at 2.27% and 59.54 keV at 35.9% branching ratio). The 
two gammas only accompany the main alpha peak (c), leading 
to alpha-gamma pileup that contributes more counts to peaks d 
and e. Values of energy and branching ratio taken from IU6I1 . 



Peak Centroid (keV) B.R.(%) 



a 


5388.0 


1.66 


b 


5442.8 


13.1 


c 


5485.56 


84.8 


d 


5511.5 


0.225 


e 


5544.5 


0.37 



(a) Top-down view 
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(b) Side view 



Figure 3: The plate that sits just atop the HPGe crystal has 
collimation holes, drilled at 0°, 30°, 45°, and 60° with respect 
to the normal of the plate, and all situated at a radius of 0.6" 
from the center of the plate (a). The alpha source is constrained 
to rotate above this circle, allowing it to shine through any one 
of these holes via an external rotational feedthrough built into 
the detector cryostat. The side view (b) shows the source and 
alpha shine through a representative hole onto the surface of the 
HPGe detector. 



at the top of the modified cryostat enclosure (b). The 
feedthrough connects to an arm and source holder (c) in which 
a windowless-alpha source is placed. The source can be rotated 
in a circle above the aluminum cap (d), through which holes 
were drilled at different angles (Fig. [3}. These angled holes 
allowed alphas of different incidence angles (0°, 30°, 45°, and 
60° with respect to the normal of the surface) to reach the crys- 
tal's (e) face. The source used was a windowless 241 Am alpha 
source from Isotope Products with an activity of 161.9 Bq [ 17]. 
The decay of 241 Am^ 237 Np (Q-value: 5637.82 keV) results in 
the emission of an alpha particle 100% of the time. The five al- 
pha branches with probabilities greater than 0.1% are noted in 
Table Q] There are two gamma emissions with non-negligible 
branching ratios that occur with the 5485.56 keV alpha. Given 
the proximity of the source to the detector, a large fraction of 
these decays will deposit both an alpha and a gamma in co- 
incidence, and the energy spectrum as measured by the HPGe 
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Figure 4: Representative energy spectrum (from simulation) 
of a collimated 241 Am source. The 5 peaks, labeled a-e, corre- 
spond to the peaks in Table Q] The centroids of the peaks are 
lower in energy than the centroids from Table Q] due to energy 
loss in the dead region of the detector, as described in the text. 



detector will reflect this with some of the events from the 5485 
keV alpha going into a higher energy peak. This structure can 
be seen in Fig. [4] 

Data were collected using four different collimation holes, 
corresponding to alphas striking the surface of the detector at 
incidence angles as described above. Figure[5]shows calibrated 
energy spectra for all four data sets. 

A charged-ion interaction model was constructed to char- 
acterize the detector's response to surface alphas. An alpha 
traversing a non-active region of length Ax will undergo a mean 
energy loss of 
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(1) 



where the stopping power (itself a function of the alpha's en- 
ergy) is integrated over Ax. A Bohr model of non-relativistic 
heavy particles on thick absorbers iflU [l9ll was used to relate 
the energy straggling to the expected width of a mono-energetic 
alpha beam incident on an absorber (i.e. a non-active region). 
An idealized energy spectrum (as measured by an HPGe detec- 
tor) of alphas of initial kinetic energy Eg traveling through an 
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absorber of thickness Ax would then be described by a gaussian 
with mean Eq - A.E and width <x: 



a 2 = AnN a r 2 e {m e c 2 fp—Ax. 



(2) 



Here, N a is Avogadro's number, r e is the classical radius of 
the electron, m e the mass of the electron, and p, Z, and A the 
density, atomic number, and atomic mass of the absorber ma- 
terial. In practice there are two further required modifications. 
Nuclear quenching within the HPGe detector results in energy 
losses that do not register with the ionization detector. This re- 
sults in a further energy offset, AE^q, that must be included. 
Furthermore, HPGe detectors do not have 100% efficiency for 
charge collection, resulting in a low-energy tail. For this reason, 
an exponentialy-modified Gaussian (a Gaussian convolved with 
an exponential) is used in the model to account for this asym- 
metry in the peak signal. Incorporating these two additions, the 
model-predicted energy spectrum from a mono-energetic alpha 
source traversing an absorber becomes 
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where //q = A£ + AE^q, t is the exponential parameter, and <x 
is defined as in Eq. [2] 

Tailoring Eq. [3] to the 24 'Am source in the test stand re- 
quires summing 5 such peaks and also taking into account the 
background from cosmic ray-induced events. The cosmic back- 
ground is well-described by a linear polynomial (Py) in the re- 
gion around the peak structure (5100-5600 keV). The energy 
loss and straggling parameters, AE and <x, depend on the path 
length through the dead layer (depth O) via Ax = £)/ cos 0, 
where 9 is the incidence angle of the alpha. The probability 
density function describing surface alphas from 241 Am imping- 
ing on an HPGe crystal surface at incident angle 6 in the test 
stand is then 



ME, 0, 



S ) = IT^ X c ^ £ ' M,<f) + — 5-*>, , (4) 

Af total f-f Aftotal 

where the coefficients represent the weights of the alpha 
peaks in Tableland sum to 1. 



The weights w; for a given single incidence angle are de- 
termined via a line-of-sight simulation that takes into account 
the collimation angle, size, and the position of the source. It is 
important to note that the peaks in Fig. [5] are wider than they 
would be with a perfect collimator. 

The analytical model was fit to the test stand data using a 
maximum-likelihood fit. All four data sets (with nominal inci- 
dence angles of 0°, 30°, 45°, 60°) were fit simultaneously with 
the goal of extracting the dead layer parameter D. The com- 
bined model for the data sets with the data are shown in Fig. 
[5] The value of the dead layer, as extracted from the fits, is de- 
termined to be 0.307 + 0.005/im. This is in agreement with the 
dead layer thickness as given by the data sheet from ORTEC 
(0.3 j-im, no stated uncertainty). 

Simulations of the test stand were also performed using the 
Geant4- and ROOT-based software package MaGe, with the test 
stand incorporated into the package's geometry. A dead layer 
was simulated in the crystal surface as a step function, ignor- 
ing any energy deposits within the dead region (this assump- 
tion is the same as our analytic treatment). The simulated out- 
put spectra were convolved with a detector-response function, 
i.e. a modified gaussian incorporating the same skew param- 
eter from the model fit. The width of the response function 
was determined by measuring the detector resolution from fits 
to gamma peaks in the detector (between 100 and 2600 keV), 
fitting the resolution values to a parametrized resolution func- 
tion f(E) = a J 1 + BE 120J], and extrapolating the resolution 
at 5.4 MeV from that function with the fitted values of a and B. 
Comparing the simulations to the data (Fig. [6} revealed discrep- 
ancies in both the width and offset of the peak structure. The 
simulated offset and width do not match with the data or the 
analytic model. A simplified simulation of a mono-energetic 
alpha incident on the surface of a detector was performed for 
a range of incidence angles (corresponding to a range of dead 
layer thicknesses) to determine the difference in alpha peak off- 
sets and widths between simulation and the model (Fig. [6). 
These discrepancies were then used to create a second convolu- 
tion function that was applied to the simulated data (as a func- 
tion of incidence angle). As shown in Fig. [5] the corrected 
simulated spectrum is in good agreement with the data and the 
analytic fit. It should be pointed out that this particular simula- 
tion was performed using Geant4 version 9.0, and simulations 
run with newer versions resulted in different discrepancies. 



2.2. Test Stand Operation 

The test stand was operated at Los Alamos National Labora- 
tory at 2200 meters above sea level. The high cosmic-event rate, 
coupled with the relatively weak source (161.9 Bq), forced the 
size of the collimation holes to be larger than ideal. The poorer 
collimation power resulted in a larger range of incidence angles 
for a particular hole than desired, and required a modified prob- 
ability density function composed of a weighted sum of single- 
angle p.d.f.s. The final p.d.f., taking into account the spread of 
angles around the nominal collimation angle 6, then becomes 



R(E, 0,D) = J^ mR(E, h D). 



(5) 



2.3. External-Bulk Backgrounds 

The rotational feedthrough was replaced with a blank-off for 
the bulk-type background study. A 4 cm length of thorium wire 
(99% 232 Th~) j21ll was looped 5 mm above the crystal face, al- 
lowing direct shine from the wire to the face of the crystal. The 
energy spectrum from this data set is shown in Fig. [7] The 
close proximity of the wire to the detector surface means that 
the detector is inundated with gammas from the 232 Th chain in 
addition to the alphas. Double and triple coincidence gammas 
from 208 T1 are evident in the energy spectrum up to 3710 keV 
(511 keV + 583 keV + 2615 keV). The step structure above 
3700 keV is evidence of thorium alphas with the high-energy 
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(a) 0° incidence. 



(b) 30° incidence. 
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(c) 45° incidence. (d) 60° incidence. 

Figure 5: Surface-type alphas from 241 Am at incidence angles of 0°, 30°, 45°, and 60° with respect to the surface normal. Also 
shown are the analytic-model fits (solid line) and simulated spectra from Geant4 (dashed line, corrected via the convolution from 
Fig.©. 



end of each step corresponding to the initial kinetic energy of a 
given alpha. 



A model, similar to that of the surface-type detector re- 
sponse, was constructed to describe bulk-type data on the detec- 
tor. The typical range of an alpha in solid materials of interest 
(thorium, copper, lead, gold, germanium) is on the order 10's of 
fim. Thus only the outer portion of material facing the crystal 
will contribute alpha backgrounds. An alpha of original energy 
Eq, emitted in an external-bulk material at a depth d from the 
material surface at an angle 9 with respect to the normal of that 
surface, will then lose an average energy similar to Eq. Q]where 
dE/dx is a function of the alpha's energy and distance traveled. 
Similarly, energy straggling will widen the resulting peak. The 
spectrum from an alpha of initial energy £o at a depth d emitted 
at an angle 9 would then be given as 



/ (E, E , d, 9) 



exp 



E + Ea-AE) 



2cr 2 



(6) 



cr £ = AnN a r e {m e c l Yp—Ax 



AE 



= 289.9Ax 

dx 



-dx. 



where we once again assume a Bohr model of energy straggling 
and o 2 is calculated using p, Z, and A of thorium. While / could 
be described with an exponentially-modified gaussian as for the 
surface data, such details are washed out in practice because of 
the continuum nature of the energy spectrum as / is integrated 
over 9 and d. 

The model spectrum for a particular decay is then found by 
integrating Eq. [6] over 9 and d. Values for dE/dx were taken 
from alpha stopping -power and range tables H22I1 . and this in- 
tegration was performed for all the alphas in the 232 Th decay 
chain. The corresponding energy spectra were then summed, 
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AE - a + pAx +yAx ! 




a = 15.59 + 0.01 




P - -7.64 + 0.05 




y = 0.50 ±0.03 
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(a) Deficit in peak offset position (AE = /i ana i - yU s j m ). 
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Figure 7: Data from a thorium wire source in the test stand. 
Alphas from the 232 Th decay chain produce steps in the energy 
spectrum. The specturm is shown only from 2.7 MeV to 12 
MeV to emphasize the high-energy alpha spectrum. Data below 
2700 keV are dominated by betas and gammas from the 232 Th 
chain. The peaks visible below 4 MeV are from double- and 
triple-coincidence gammas from 208 T1. 
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- a = 0.17±0.06 




~ p= 112.59±0.20 
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(b) Differences in peak width between simulation and the model. A paramet- 
ric function of dead layer distance (<r = sja + fiSx is shown for the analytic 
model (red, dashed line) and is also fitted to the simulated points (blue, solid 
line). 

Figure 6: Discrepancies between simulation and analytic 
model. When compared to the measured data, the Geant4 sim- 
ulation overestimates the peak position offset (a) while under- 
estimating the peak width (b). A quadratic function was fit to 
the peak position offset and a square root function was fit to 
the width offset. These functions were then used to generate a 
convolution to modify the simulated energy spectra. The off- 
set deficit at distance is consistent with the expected offset 
from nuclear quenching which is not taken into account in the 
original Geant4 simulation. 



assuming the chain is in secular equilibrium. This assumption 
is valid in our comparison with our data because all alpha de- 
cays in the 232 Th chain, except the first, happen within days of 
each other. The first alpha decay in the chain ( 232 Th) emits a 
4 MeV alpha which is not discernable below the beta/gamma 
continuum in the data set. 

This model was compared with the data set (Fig. [SJ. The fast 
alpha decay (300 ns) of 212 Po, coming after the beta decay of 
212 Bi, occasionally results in pileup in the detector of the alpha 
(8.8 MeV) and the beta (endpoint 2252 keV, intensity 55.4%). 
To treat this, the p.d.f. representing the 2I2 Po alpha spectrum 
was convolved with the beta spectrum of 212 Bi. A combined 
p.d.f. consisting of the unmodified 212 Po p.d.f., the 212 Po+ 212 Bi 
p.d.f., the upper 232 Th chain p.d.f., and a cosmic-ray spectrum 
p.d.f. consisting of a constant plus exponential function was 
constructed. This combined p.d.f. was fit to the data set and the 
results are shown in Fig. [8] 



8000 



u 6000 
M 
o 

<r> 

2 4000 
c 

o 
U 

2000 









Bulk-Alpha Test Stand Data 




Alpha Model + Cosmic BG 




Upper 232 Th chain a s 




212 Poa 




2I2 Poa • Hi |1 




Cosmic BG 







4000 5000 6000 



7000 8000 9000 
Energy [keV] 



10000 11000 12000 



Figure 8: Fit of bulk-alpha energy spectrum with composite 
model from a thorium-wire source in the test stand. 
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3. In Situ Measurement 

The model developed in the previous section was used to an- 
alyze the energy spectrum of in-situ backgrounds from an ex- 
isting underground detector. An n-type HPGe detector, under- 
ground at the Waste Isolation Pilot Plant near Carlsbad, NM, 
was discovered to have a peak in its background energy spec- 
trum at 5.3 MeV as seen in Fig. [9]. We model the spectrum 
assuming that all events above 2.7 MeV (i.e. above the 2615 
keV gamma peak from the decay of 208 T1) must originate from 
either surface-alpha events, external-bulk alpha events, or cos- 
mic events. A muon-veto scintillator panel was installed above 
the detector for a portion of the data taking, giving an energy 
spectrum of cosmic-ray muons from events that were tagged 
with the veto. The same constant plus exponential function 
that was used for the test stand data in Sec. [2] was fit to this 
cosmic-tagged data and subsequently used in the background 
fit. A p.d.f. — consisting of individual p.d.f.s of the surface- 
and external-bulk 232 Th, 238 U, and 210 Po events and the cos- 
mic p.d.f. — was used to fit to the spectrum. This fit is shown in 
Fig. [10]and the results are summarized in Table|2] Components 
from the 238 U chain — save 210 Po — were consistent with zero 
and so were excluded from the final fit. Alpha yields and rates 
are tabulated in Table [2] The only evidence of bulk events from 
the 238 U chain were found to be 210 Po. This could conceivably 
come from bulk material containing short half-life daughters 
lower in the 238 U decay chain (anything below 226 Ra). This 
would have decayed to the long-lived 21() Pb (22 year half-life), 
allowing a lead-supported supply of 210 Po alphas. 
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Figure 9: Energy spectrum from underground HPGe detector 
showing a 5.3 MeV alpha peak (characteristic of the alpha from 
210 Po decay). 



4. Alpha Backgrounds and Ov/?/? 

The efficiency for a surface alpha decay on a p + surface, 
e.g. the alpha from 210 Po, to populate the region-of-interest 
in a double-beta decay experiment using HPGe detectors was 
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Figure 10: Fit to the spectrum of an underground HPGe de- 
tector's energy spectrum with a composite alpha model. The 
largest contributions above 2700 keV are from 210 Po surface 
events and from external-bulk 232 Th and 210 Po events. 

Table 2: Yields and rates of individual components of com- 
posite alpha model, fitted to high-energy data of underground 
HPGe detector. The data set represents 3 17.27 days of livetime. 



Contribution 



Yield Counts / Day 



Cosmic Events 


1033+g 


T TO+0.21 
J.JJ_ QA9 


Surface 210 Po 






Bulk 210 Po 


641$ 


1 89 +0 - 31 


Bulk 232 Th 




a 97+O.6I 



calculated from simulations using MaGe. For a nominal dead 
layer value of 0.3 //m, the probability that the energy detected 
from a 2I0 Po decay on the surface falls within the Qv/3/3 region- 
of-interest (ROI) is (1.21 ± 0.05) x 10~ 5 . Assuming a possi- 
ble variation in dead-layer of +0.01 fim adds a systematic un- 
certainty of ±0.04 x 10~ 5 . This includes the corrective con- 
volution that was required in Section 12.11 although the differ- 
ence in calculated efficiency with and without the convolution 
was only 0.01 x 10~ 5 . The final, simulated efficiency is then 
(1.21 ± 0.05(stat)±0.04(sys)) x 10~ 5 . The same efficiency was 
also calculated using the surface alpha analytic model. The 
model was used to generate an energy spectrum from surface 
210 Po decays, and the efficiency for a decay to land in the 4 
keV-wide ROI is calculated to be (1.467 ± 0.004) x 10~ 5 . Al- 
lowing for the same ±0.01 /mi variation in dead layer, the model 
efficiency becomes (1.47 ± 0.004(stat)+° ^(sys)) x 10~ 5 . 

The two techniques differ by 0.26 x 10~ 5 . The test-stand data 
was in better agreement with the analytic model — particularly 
the bulk comparison — and so this value is used as the final 
calculated efficiency. The discrepancy between simulation and 
model is added as a systematic lower uncertainty. Another pos- 
sible systematic effect arises from uncertainty in the dead-layer 
profile. Both the simulation and the analytic model assume a 
step-like efficiency function for the dead layer. In this case, 
no charge collection happens at all within the dead region, but 
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charge collection is 100% efficient within the active region. The 
stopping power, 4|, is a non-linear function of alpha energy, 
and so the amount of energy lost in a dead region is depen- 
dent upon both energy of the alpha and the charge-collection 
efficiency within that region. A smoothly-varying dead-layer 
profile will result in less charge being collected than a step-like 
profile. Assuming a worst-case scenario (with the efficiency 
profile a linear function, representing the largest possible dif- 
ference in the energy spectrum), the energy spectrum can be 
compared with the uncorrected (step-function efficiency pro- 
file) and this is shown in Fig. QT| When this correction is added 
to the simulated and model energy spectra, the efficiency for a 
decay of 210 Po to land in the ROI increases by 7%. A related 
effect would come from non-uniformity in the dead layer over 
the surface of a detector. If the dead layer varies from point 
to point on the detector, then the efficiency for an alpha decay 
would vary from point to point as well. Accounting for the dif- 
ferent predictions of the simulation vs. the analytic model, and 
folding in the possibility of a dead-layer effect, the value for the 
efficiency of a 210 Po alpha, emitted from the p + (thin) surface 
of an HPGe detector, is (1.47+°-*° x 10~ 5 ). 

The goals of the Majorana Demonstrator will be to test the 
OvjGyS claim by Klapdor et al. [23] and to demonstrate the back- 
ground goals necessary for a future tonne-scale experiment. 
The alpha background rate of a detector (or array of detectors), 
in units of background counts within the 0v/3/3 ROI per tonne- 
year, is given as 



R 



a = ^ f 3l(r)s(r)Cl(r)dS , 
M Js 



(7) 



where M is the mass of germanium in the detectofl J{ the 
surface-alpha activity in Bq/cm 2 , and s(r) the efficiency for a 
decay at position r to count as a background. The integral is 
over all relevant surfaces dS , and Q(r) is the solid angle of p + 
area with which the area element dS has direct line-of-sight. 
The coefficient k converts the subsequent background rate from 
counts per second-kilogram into counts per tonne-year. This 
formula can be significantly simplified for decays that occur on 
the surface of the detector itself: 
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Figure 1 1 : The dead-layer profile (as a function of depth) af- 
fects the energy spectrum for surface alphas emitted from 21l) Po 
(5.3 MeV). The efficiency for surface decays is shown over (a) 
the full simulated energy range and (b) a close up around the 
region-of-interest at 2039 keV. The solid, black line shows the 
energy spectrum from the analytic model with a step-like dead- 
layer profile. The red, dashed line assumes a linear profile. 
This profile represents the greatest difference in energy spec- 
tra, and the resultant change in efficiency for a surface-alpha 
background hit is 7%. 



Ra — k j^^avgE — k\$K aV gt 



(8) 



with S the total susceptible surface area, J{ avg the average 
surface-activity rate, and e the efficiency calculated in the pre- 
vious section (a true "surface" event with solid angle Q = 2n). 
The factor A is the ratio of susceptible surface area to active 
mass (with units of area/mass). Because the Demonstrator will 
be composed of modified, p-type point contact (p-pc) detec- 
tors, it will have a particularly-low susceptibility to alpha back- 
grounds, especially in comparison to n-type detectors (Table|3]l. 
This is because the vast majority of surface area that is n + , or 
"thick". The susceptibility for the modified p-pc detectors ex- 
pected to be used in the Demonstrator is 0.34±0.03 cm 2 /kg, 



'The usual method is to use M as the amount of active mass, but alpha 
backgrounds will not differ for enriched vs. unenriched detectors. 



where the uncertainty comes from the quoted tolerances of the 
detector dimensions for the Canberra BEGe " fl 'Ge detectors that 
will be installed in the first module of the Majorana Demon- 
strator. For comparison, the susceptibility for a typical n-type 
detector is over 400 times greater. Plugging in the susceptibility 
and the calculated efficiency, the rate then becomes 



= 1.571° # X lOXvg, 

1+0.20 , 
--0.21 1 - 



J^avg in Bq/cm 
J{ avg in Decays/Day /cm 2 (9) 



Table |4] displays expected background rates from 210 Po alphas 
for several values of surface activity 3l avg . The alpha rate from 
the Sudbury Neutrino Observatory's 3 He proportaional counter 
detectors represents the cleanest detector surface ever measured 
in terms of alpha contaminationj^- Also shown are the rates 
for a typical n-type detector, given the same activity rates. 
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Table 3: Surface areas categorized by surface type for a typical n- and p-type detector, and the modified p-pc natural germanium 
detectors used in the Majorana Demonstrator. The factor A is the alpha-susceptibility factory as defined in the text. 



Detector Type 


Mass [kg] 


Surface Area/Ratio [cm 2 /%] 




A [cm 2 /kg] 






passivated p + n + 


Total 




n-type 
p-type 
p-pc 


1.1 
1.1 

0.579 


30.4/12.7 177.5/74.1 31.8/13.3 
1.8/ 0.8 35.4/ 14.8 202.6/84.5 
1.9/ 1.8 0.2/ 0.2 104.5/98.0 


239.7 
239.8 
106.6 


161.4 
32.2 
0.35 



If the alpha decays are not coming from the surface of the 
HPGe crystal, then the solid angle simplification no longer ap- 
plies. For configurations involving complicated surfaces facing 
the p + area of a detector, the integral in Eq. |7]can be calculated 
using Monte Carlo methods. An upper limit on the activity can 
also be placed, assuming that the p + region of the crystal has a 
line-of-sight view of a contaminated surface with average sur- 
face activity J[ aV g,ext- The limit is then 



Ra,ext — k^^avs.extE 



(10) 



where the equality holds if all of the solid angle that the sen- 
sitive area "sees" is emitting alphas at the surface rate 3\. a vg,ext 
(and Eq. [8]is the limiting case for this). Equation [TOl holds via 
a simple flux argument. 



Table 4: Estimated background count rates from surface alphas 
for modified p-pc detectors as will be used in the Majorana 
Demonstrator, and for typical n-type detectors. These rates 
are based upon assuming various surface-activity levels. The 
clean room entry corresponds to a class- 100 clean room built 
for radon-deposition testing for the Borexino experiment [13], 
and surface-activity values in the table represent the radon- 
deposited activity on nylon in that clean room. The calculation 
for the SNO 3 He detectors is found in 11411 . 



Source 


Activity 




Background Rate 




[Bq/cm 2 




[Counts in ROI / t-y] 








P-PC 


n-type 


Clean Room 


1.0 x 10 


-6 


0.16 


70 


In Situ Detector 


9.0 x 10 


-7 


0.14 


61 


MJ BG Model 


5.0 x 10 


-7 


0.08 


35 


SNO 3 He Detector 


5.0 x 10 


-9 


0.0008 


0.35 



Majorana from a surface-alpha standpoint (to say nothing of 
their other important qualities). As Table|4]notes, the difference 
in surface-alpha background rates between p-pc detectors and 
n-types is large. Without further cuts from pulse-shape analy- 
sis, Table [4] makes it clear that n-type detectors are unsuitable 
for Qv/3/3 experiments without heroic measures to limit surface 
activity. A final comment should be made about the passivated 
surface-the annulus of material between the n + and p + surfaces. 
This area has 10 times the surface of the p + layer for the mod- 
ified p-pc detectors, and so must be addressed. Our test stand 
was modified to allow the 241 Am source to have a direct shine 
path on the passivated surface of the crystal. Not only were no 
alphas present visible in the resultant energy spectrum, but no 
59 keV gammas were visible either. If the passivated surface 
were as susceptible to alphas as the p + surface, the background 
count rate from a given surface contamination would increase 
by a factor of 10, but still be far better than that from n-type 
detectors. We are continuing to study the issue. 



5. Conclusion 

We built a test stand for the purpose of studying surface-alpha 
decays of radioactive isotopes on surfaces of HPGe detectors. 
A simple detector response model was constructed and fit to 
test-stand data, and Geant4 simulations were compared with 
test stand data and the detector response model. The efficacy 
of simulations of surface-alpha decays for Ov/3/3 experiments 
was examined and we propose a correction factor to be used for 
such simulations. We calculate allowable contamination values 
for n-type and p-pc HPGe detectors used in double-beta decay 
searches given desired surface-alpha background assumptions. 



The susceptible surface of a modified p-pc detector, located 
right at the point contact, only has a direct line-of-sight with 
the detector mount. Alphas can only pose a background if they 
originate from surface plate-out (on the surface of the detector 
or the detector mount) or from the bulk material of the detector 
mount. Because of this compartmentalization, the rate of alpha 
backgrounds (in counts per tonne-year) should be the same re- 
gardless of the number of detectors. The rate formula for p-pcs 
(Eq. [9} would still be valid, then, for a one-tonne scale exper- 
iment made up of p-pc detectors with the same susceptibility 
factor A. The usage of p-pc detectors is extremely beneficial to 
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